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INTRODUCTION 


A radiometer is a system designed to measure the electromagnetic radiation re- 
sulting from thermally-induced molecular motion in matter. Radiometers are basically 
power-measuring receivers in which the most critical performance parameters are gen- 
erally held to be sensitivity and linearity. Man-made radio-frequency interference 
(RFI) constitutes a potential error source for microwave radiometry that is generally 
underestimated. RFI is an insidious problem in that it may be impossible to know 
that interference is present until it becomes so large that the radiometer output is 
grossly erroneous. 

RFI may be classified as being either in-band or out-of-band relative to the 
radiometric channel bandpass. While this paper is concerned with minimizing the 
deleterious effects of out-of-band RFI on radiometer performance, it is worth noting 
that it may be possible to establish the presence of RFI, of either category, by 
observing any characteristics of the interference which differ from the geophysical 
signal. In this manner, it should be possible to achieve some degree of data 
correction. 

Out-of-band interference can produce radiometric errors through several system 
imperfections. The first and most obvious is non-ideal, i.e., not rectangular, 
channel selectivity wherein one or more nearby but still out-of-band interfering 
signals reach the detector along with the desired signal. This is a linear effect 
which can, in principle, be reduced to any desired degree by improved selectivity. 

In practice, sensitivity considerations usually result in the bulk of the selectivity 
being located well down the signal processing chain (e.g., after the amplifier and 
miKer stages) where bandpass filter insertion loss does not significantly degrade 
system sensitivity. However, this arrangement leaves the system elements preceding 
the primary selectivity susceptible to the nonlinear effects of desensitization and 
in-band spurious generation from out-of-band RFI. These nonlinear effects are ir- 
reversable unless the RFI is clearly identifiable so that it can be "removed" from 



the data. Prevention is much preferable to after-the-fact data correction. However, 
prevention by preceding potentially nonlinear front-end elements, such as amplifiers 
and mixers, with additional selectivity incurs a cost in sensitivity. 

This study examines the trade-offs between sensitivity and interference immunity 
associated with front-end bandpass filtering. The results are presented in a form 
that defines the radiometric sensitivity degradation factor, AT'/AT, in terms of 
various filter characteristics. Additionally, the "optimum" filter configuration can 
be totally defined once the independent (filter) variables of percentage bandwidth, 
q-factor(s), and desired band-edge rejection are specified. 

LIST OF SYMBOLS 

Aq numerator coefficient of t(s) 

Ap numerator coefficient of s ) when predistortion is applied 

a magnitude of real part of poles closest to jar-axis (in main body) 

a attenuation constant (in Appendix A) 

B bandwidth of bandpass filter between attenuation points 

B no radiometer output low-pass noise bandwidth 

B r arbitrary reference bandwidth 

B statistical bandwidth 

b 

B g £ statistical bandwidth of predetection input filter 

0 phase constant 

C constant modifying gain of predistorted filter realization (C < 1) 

d dissipation factor 

AT unmodified radiometer sensitivity 

AT' modified radiometer sensitivity with degradation 

f frequency variable, Hz 

f Q band-pass center frequency 

fj lower band edge frequency 
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upper band edge frequency 
radiometer amplifier power gain 
net input reflection function 
load reflection coefficient 
net output reflection function 
source reflection coefficient 
RMS load current 

Boltzmann's constant = 1.38 x lO -2 '* J/K 
power dissipative loss function 
integrated loss factor 
frequency variable, rad/sec 
imaginary part of an arbitrary pole p^ 
available source power 
an arbitrary pole of t(s) 
filter output power 

low-pass unload element quality factor 
band-pass loaded resonator quality factor 
passband ripple in dB 
bandedge attenuation in decibels 
load resistance 

image impedance when it is real 
source resistance 
Laplace variable 

real part of an arbitrary pole p^ (ctr < 0) 
available power density of source 
incident power density at filter input 
output load power density 
available internal noise power density 



t(s) 

t'(s) 


t d (s) 

t^( j 0)) 
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a INJ 

T ' 
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l REF 
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z ln ( jw) 
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transmission function 
modified transmission function 
transmission function with dissipation 
predistorted lossless transmission function 
integrated transmission factor 
antenna brightness temperature 
injected noise 

modified injected noise after filter is added 

filter output noise temperature 

physical temperature of filter 

receiver noise temperature 

reference noise temperature 

equivalent source noise temperature 

system noise temperature 

image propagation constant 

RMS source voltage 

short-circuit filter transfer admittance 
input impedance of filter 
image impedance 
output impedance of filter 


NOISE MODEL FOR FILTER 

The manner in which the noise introduced by the dissipative elements within a 
filter interacts with the radiometer input signal will be developed in this section. 
Consider the block diagram of a filter with source and load terminations as shown in 
figure 1. Assume that all dissipative elements within the filter are at a physical 
temperature T p . 
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The frequency-dependent reflection functions at the input and output are denoted 


as Yi(jw) an d y 0 C j cu) respectively. These functions represent the net effect of 
all multiple reflections and could be expressed in terms of the simpler frequency 
independent terminal reflection coefficients based on the source and load resistances 
and the image impedance. However, it is not necessary to express an image impedance 
explicitly, and this is rarely done with s-plane filter designs, such as were used in 
this study. A development of the relationship between the terminal reflection coef- 
ficients and these net functions is given in Appendix A. 

Let t(joj) denote the transmission function defined on a magnitude-squared 
basis as: 


tfi )|2 output load power density 
J ' available source power density 


( 1 ) 


In the development that follows, it is convenient to define a modified 

A 

transmission function |t'(jw)| as 


t' (ja))| 2 


output load power density 
incident power density 


( 2 ) 


where the incident power density is the fraction of the available power density 
accepted by the filter. 

Let S av (f) represent the source available power density, and assume that it is 
flat over the frequency range of interest and therefore representable in terms of an 
effective source temperature T g . Thus, 

S (f) = kT W/Hz (3) 

av s 

where k = Boltzmann's constant = 1.38 x 10"^ J/K. 


Let S^(f) and S Q (f) represent the incident and output spectral densities 

2 

respectively. The power density reflected at the filter input is | y. ( jcu) | S 

1 3V • 
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Thus, the incident power spectral density is related to the available power density 
by 


S ± (f) " [1 - lYiCju,) | 2 ]S av (f) 


(4) 


The modified transmission function is 


1 1 ' ( joj) | 2 


s o ( f ) 

W [1 


S o( f > 

|Y 1 (ja> ) | 2 ] S av (f) 


(5) 


The transmission function is 




s o (f) 

s „< f > 


( 6 ) 


Comparing (5) and (6), the following relationships between |t'(ju>)| 2 and |t(ju))| 2 
may be deduced: 


t ' ( jm) | 2 


1 1( jto) 1 2 
1 - I y ± C ju>) | 2 


(7) 


or 


| t(j o>) | 2 = [1 - |y i (jo))| 2 ] |t’(jio)| 2 (8) 

The manner in which noise associated with the dissipative elements within the 
filter interacts with the input noise is illustrated by the block diagram of 
figure 2. Let Sp(f) represent the available internal noise power density, and 
assume that 

S (f) = kT W/Hz (9) 

P P 
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The noise power propagates in both directions. At the input, the source resistance 
accepts from the filter a power density [l - | Y ^ ( J w ) | ]S ^ ( f ) - For thermal equi- 
librium, this power must be subtracted from the source power, so the net power 
accepted by the filter from the source is 

[1 - |y ± ( jo>) | 2 ]S av (f) - [1 - | Y 1 ( jw) | 2 ]S p ( f ) 

as shown to the right of the input summing junction. This power density is altered 
by the modified transmission function. The effect of the output mismatch on this 
noise has already been considered in the definition of Yj(joi). 

In addition to internal noise propagating to the left as already considered, 
internal noise propagates to the right. At the load junction, power density 
I Y 0 C j to) | ^ S (f) is reflected back to the network, so the net contribution of the 
internal noise to the load is [l - | y q ( jw) | 2 ]Sp(f) . 

Combining all the preceding sources, the net output spectral density S Q (f) can 
be expressed as 

S o (f) = {[1 - | Y± ( j^o) | 2 ]S av (f) - [1 - |Y i (jo))| 2 ]S p (f)}|t*(j w )| 2 

+ [1 " lY 0 (ju>)| 2 ]S p (f) (10) 

As shown underneath the block diagram in figure 2, this expression can be simplified 
to 

S Q (f) = 1 1( jco) | 2 S ay (f) + [1 - 1 1( j o>) | 2 - | Y 0 ( j to ) | 2 ] S p (f ) 

A loss function |£(ju))| 2 will now be defined as 

| £( joi) | 2 = 1 - 1 1( jtu) | 2 - | Y q C j to ) | 2 
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(ID 


( 12 ) 



It is interesting to note that while the focus of analysis is on transmission from 
source to load, the pertinent reflection function is the output function. However, 
for networks initially designed on a lossless basis with resistive termination at 
both ports, it is shown in Appendix B that y 0 ( j to) = YfOw), so either port could be 
used to establish the loss function in that case. This situation is not true in 
general. 

Utilizing the preceding definition, the output power density can be expressed as 


S Q (f) = |t(jo))| 2 S av (f) + |£(ja))j 2 S (f) (13) 

INTEGRATED FILTER FUNCTIONS 

Based on the preceding development, assume that the differential output power 
AP q in a differential bandwidth Af is to be determined. This quantity is readily 
expressed as 


AP„ = 1 1( jco) | 2 S (f ) Af + | £(joj)| 2 S (f)Af 

U d V p 


= kT j t( jrn) | 2 Af + kT U(jo))| 2 Af 
s p 


(14) 


The total power P Q in a bandwidth B = f 2 - fj is determined by changing the 
differential formulation to integral form, and the result is 


P 

o 




2 

U(jw)| df 


(15) 


It is convenient to express the result in terms of an effective temperature T Q , 
defined over a rectangular bandwidth B according to 

P = kT B (16) 

o o 
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Substitution of (16) in (15) and cancellation of common factors yield 


T 

o 


if 2 1 1( jio) i 2 df 


+ T 


i f 2 U(ja))| 2 df 


( 17 ) 


2 2 

The integrated transmission factor x and the integrated loss factor L will 

O 

be defined as the symbol t in the definition of t should not be confused with 
the integration time parameter appearing in radiometer sensitivity equations. 


t = 


if 2 |t(jm)| 2 df 


i 2 = i f 2 ua<o)i 2 df 


(18) 


(19) 


With these definitions, T Q can be expressed as 

T = t 2 T + L 2 T (20) 

o s p 


DEGRADATION FACTOR 

The preceding definitions will now be used to determine the degree to which the 
interference-rejection filter degrades the radiometer sensitivity. This development 
is limited to the closed-loop noise injection feedback radiometer, the configuration 
of primary interest. The corresponding effect in a total power radiometer is sub- 
stantially different. 

The interference rejection filter is located after the directional coupler and 
is designed to provide a specified attenuation R (in dB) at two specified band- 
edges. The basic form of a typical response is shown in figure 3. It is assumed 
that the primary predetection bandwidth selectivity is obtained after several stages 
of gain and has a much more ideal, i.e., rectangular, response over B = f 2 - f j . 
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The attenuation of the interference rejection filter will thus be specified at these 


frequencies, and the amplitude response will be assumed to be truncated abruptly as 
shown by the dashed lines. This realistic assumption considerably simplifies the 
integration process as will be demonstrated later. 

The pertinent portion of the system affected by the filter is shown in figure 4. 
The system without the filter is shown in figure 4(a). Sufficient injected noise 
temperature T IP jj is added to the antenna temperature T A such that 
T A + T IN j = T ref , where Tre F is the reference temperature. This power balancing 
condition is maintained by the closed loop. 

When the dissipative filter is added as shown in figure 4(b), the resulting 
integrated transmission function is constrained by the inequality t < 1. The 
injection noise Tj.^j increases to a level such that the filter output equals 
Tref as before. Thus, the effective overall noise level at the Dicke switch is not 
changed. However, the loop gain has been reduced by t . In order to maintain the 
same loop time response and low-pass output noise bandwidth, a gain of 1/x must be 
added to the loop as shown. Since 1/x^ > 1, the fluctuations within the loop are 
accentuated by this process, and the sensitivity is degraded. 

The second effect of the filter is to reduce the statistical bandwidth. The 
statistical bandwidth B g in general is defined by 


B 

s 



|t(ju))| 2 df 
|t(ja))| 4 df 


( 21 ) 


This quantity is the bandwidth that should be used in the radiometer sensitivity 
equation. Based on the assumed ideal overall predection filter characteristic, the 
statistical bandwidth used in this study is therefore: 
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r / 2 

|t(ja))| 2 df 

L f l 


r f 2 
/ 2 

|t(ja>)|* df 

L f i 



( 22 ) 


Let AT represent the idealized sensitivity of the radiometer without regard to 
the interference filter. This widely used definition is 


41 ■ 2i sys 



(23) 


where the system temperature T SYS is T SYS “ T R + T REF and T^ is the receiver 
temperature. The quantities B^ and B no represent, respectively, the detector 
input statistical bandwidth and the output noise bandwidth. For the ideal rect- 
angular bandpass B^ = B. Thus, without the interference filter, the sensitivity is 


AT 


2 T 


SYS 



(24) 


Let AT' represent the sensitivity as degraded by the interference filter. 

Since the effective statistical bandwidth is now B g , and the noise fluctuations are 

2 

accentuated by a factor 1 /t , equation (23) is modified by 


AT' 


2 


T SYS j 2B no 
t 2 B s 


Division of (25) by (24) results in 


AT' 

AT 



(25) 


(26) 
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This ratio will be referred to in all subsequent work as the degradation 


factor . It is a measure of the degradation in the sensitivity resulting from the 
loss in gain and the non-rectangular bandshape response of the interference filter, 

FILTER FUNCTIONS 

The actual interference rejection filter used in a radiometer is normally a 
band-pass filter. However, the vast majority of band-pass filters are designed from 
low-pass prototypes using the low-pass to band-pass geometric transformation. 

Further, the analytical properties of the band-pass response may be mapped to the 
low-pass region with a resulting one-to-one correspondence between the regions. 

The approach used in this study was to perform all the analysis with the low- 
pass prototype functions. As previously mentioned, the various properties tabulated 
are invariant with respect to the transformation. The number of poles listed for the 
low-pass prototypes must be interpreted as the number of resonat o rs (or, equiv- 
alently, the number of pole-pairs) when translated to the band-pass case. 

The number of low-pass poles was varied from two to seven in the study, thus 
allowing band-pass filters with from two to seven resonators to be included. Filter 
types included Butterworth and Chebyshev with ripple ranging from 0.1 dB to 3 dB. 

The effects of losses were included, and both non-distorted and pre-distorted designs 
were considered. 

Let t(s) represent the s-domain transmission function. For a low-pass filter 
with all transmission zeros at s = °°, this function for n poles may be expressed 
as 


t(s) 


(s-p 1 )(s-p 2 ) 


(s-p n ) 


(27) 


12 



where A Q is a constant and the various Pk = CT k + represent the poles. The 
quantities are the real parts (a^ < 0) and are the imaginary parts. The 
poles and constants of the filter types considered are tabulated in Appendix C. 

Assume that there is uniform dissipation occurring in all the band-pass 
resonators. In the low-pass prototype, this is accounted for by defining a 
transmission function t s ) with dissipation as 

t d (s) = t ( s + d) (28) 


where d is an appropriately defined dissipation factor. With non-predistorted 
filter designs, the dissipation factor is applied as given in (28). 

When predistortion is desired, a modified lossless transmission function 
t^(s) is first defined as 

t^(s) = t ( s - d) (29) 

which amounts to shifting all poles of t(s) to the right. The maximum value of d 
is limited to the smallest real part of any of the poles. This function is then 
realized as a lossless network. Because of maximum power constraints, the practical 
realization of the associated peaked response results in Ct(s - d), where C < 1. 
When dissipation is subsequently applied, the function ( s ) now becomes 

t,(s) = Ct.(s + d) = Ct(s + d - d) = Ct(s) (30) 

d x 


The effect, the desired passband response shape is preserved through the process. 
Unfortunately, predistorted designs suffer more midband loss as a result of the pole 
shifting, with possible deleterious consequences in radiometer applications. 

The magnitude-squared transmission function for any t(s) may be expressed as 


| t(jco)| 2 = t(s) t(-s) 


(31) 
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o 

The resulting function is a rational function of w and may be easily formulated on 
a geometrical basis in terms of the pole locations, the frequency w, and the dis- 
sipation factor d. 

For predistorted designs, a significant loss occurs because the initial shift of 
the poles to the right creates a more underdamped response with significant peaking. 
Let A q represent the numerator constant of t(s). It can be shown that a very good 
approximation to the corresponding constant Ap for the predistorted filter is 

A = (— — d ) A (32) 

p V a / o 

where a is the magnitude of the real part of the pair of poles closest to the 
joj-axis. 


COMPUTATION OF DEGRADATION FACTORS 

Several computer programs were developed for evaluating the degradation factor 
as a function of the number of poles, the pole values, dissipation losses, and band- 
edge attenuation. The various programs evolved over a period of time as data were 
plotted and evaluated. 

The two major programs used in the final tabulation of data are defined as 
"FILTER1 " and "FILTER2", and they are written in BASIC. Listings of these programs 
are provided in Appendix D. 

Both programs perform a numerical integration of the transmission function 
magnitude-squared and magnitude to the fourth power so that statistical bandwidth and 
degradation factor may be computed. 

PRESENTATION OF GENERAL RESPONSE CURVES 

The program "FILTER1" was used to generate a large quantity of data for various 
filter types, order, attenuation, and dissipation factor d. The results are shown 
graphically in Appendix E. Each figure represents a different filter order and type. 
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Let q represent the unloaded "Q-factor" for the normalized low-pass prototype 
design. Since w = 1 rad/s was assumed as the normalized reference, 



(33) 


The dissipation factor was used as the parameter in generating the curves, but 
was converted to the equivalent q value for plotting. 

The band-pass filter Q is related to the low-pass prototype q by 


f 



where f is the band-pass center frequency and is the low-pass reference 

bandwidth. 

Although the data of Appendix E are useful and display some interesting trends, 
there is a subtle point that must be considered in interpreting the results. The 
value of d (and subsequently of q) is based on the bandwidth of the reference 
definition of the particular filter type. For example, assume that a 0.5-dB 
Chebyshev filter is to be modified with a dissipation d = 0.1. This means that 
q = 10 at a frequency to = 1 rad/s in the low-pass prototype. Unless the design is 
predistorted, however, there is no simple way of predicting the exact shape of the 
response after dissipation is added. This poses some difficulty in making a "fine" 
interpretation of the results of Appendix E. 

A somewhat different approach was used in the development of the program 
"FILTER2", and these results are tabulated in Appendix F. All of the curves obtained 
in these presentations were based on a specified attenuation of either 10 dB or 
20 dB. The two frequencies producing these values of attenuation with dissipation 
added are determined numerically. The value of q is adjusted so that the value 
plotted corresponds to the exact q at the bandedge where the desired attenuation is 
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achieved. For example, with R = 20 dB and q = 30, the low-pass prototype filter 
has q = 30 at the exact frequency where R = 20 dB. 

For Appendix F, therefore, the relationship between low-pass prototype q and 
band-pass Q is 

f 

Q - q (35) 

where B is either the 10-dB or 20-dB bandwidth. 

Some summary results are provided in the three figures of Appendix G. Figures 
G1 and G2 provide some curves of degradation factor versus the number of poles with 
q as a parameter for 10 dB and 20 dB respectively. Although the independent vari- 
able is obviously restricted to integer values, the curves are enhanced by 
extrapolating between points. 

Figure G3 provides a macroscopic range of operating conditions between three 
poles and seven poles with low-pass q between 30 and «>. The lowest order (two 
poles) was not used because it is felt that three poles would represent the lowest 
order of practical interest. 
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SUMMARY AND CONCLUSIONS 


This paper quantifies the sensitivity loss which results from the presence of 
interference-rejection bandpass filters in the front-end of noise injection feedback 
radiometers. This sensitivity loss was expressed in terms of a "radiometric sensi- 
tivity degradation factor", which is the ratio of the actual sensitivity AT' to the 
theoretical sensitivity AT. Sensitivity degradation results from the combined 
effects of changes in the input statistical bandwidth and changes in the transmission 
function resulting from the filter. Both of these changes are determined from 
integration of the appropriate frequency-dependent functions. 

In the course of the investigation, various models concerning noise and dissi- 
pative effects were developed. Some of these models show promise for additional 
investigations on radiometer sensitivity. Computer programs employing numerical 
techniques were developed for processing the data. Numerous data were generated and 
plotted. Several common filter types were employed and various unloaded resonator 
quality factors were assumed. This resulted in a comprehensive compilation of data 
having potential use in radiometer system design and analysis. The results of this 
paper allow the radiometer system designer (a) to make an intelligent choice with 
regard to the number of poles, and passband ripple in the front-end filter, given the 
inherent Q of the filter elements, and (b) to recognize the radiometric penalty 
incurred by the interference immunity gained with the input filter. 
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[l-|y (ju))| 2 ]S (f) 

X dv 

-[l-|Y.(jco)| 2 ]S p (f) 



S o (f) = { [1-1 Y ± (jco) [ 2 ]S av (f) — 1 1— I Y ± ( jto) I 2 ]S p (f) 1 1 ’ (jto) [ 2 +[ 1- 1 Y Q ( jw) | 2 ] S p ( f ) 

= [l-|Y,(j«)| 2 ] |t'(jo))| 2 [S (f)-S (f)]+[l-| Y (joo) 1 2 ]S (f) 

X a. v P Op 

= | t ( jw) | 2 [ S ( f ) — S ( f ) ]+[ 1- | Y ( j“) i 2 ] S (f ) 

clV p Op 

= 1 1 ( j oj) | 2 S ( f ) + [ 1— | t ( jaj) | 2 -| Y (jw) | 2 ] S ( f ) 

cL V Op 

= | t(jw) | 2 S (f) + I 1( jw) I 2 s (f ) 

dv P 

Fig. 2. Block diagram of filter showing interaction 
of dissipative losses and input noise. 






Typical response of interference filter coupled 
with more ideal detector input response. 
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Fig. 4. Block diagram of portion of radiometer showing 
addition of filter and compensation. 






APPENDIX A 


INTERACTION OF REFLECTIONS 

Consider the transmission system shown in figure A1 with a source resistance 
R g and a load resistance R^. Using image parameter network theory, the image 
propagation constant is 0 = a + jg, where a is the attenuation constant and 0 
is the phase constant. The image impedance is assumed to be Z Q . The source reflec- 
tion coefficient y s is defined as 


Y 


s 



R + Z 
s o 


(A-l) 


The load reflection coefficient Yl *- s defined as 


Y t = 


R T - Z 
L o 


L R L +Z o 


(A-2) 


Let Z^ n ( j to ) represent the complex impedance reflected to the input. This 
function can be shown to be 


z in(j^ = z 0 


1 + Y L e 


-26 


1 ' 


-26 


(A-3) 


Let z 0 ut<J oj) represent the complex impendance reflected to the output. This 


function is 


Z (jo,) = Z 
out J o 


1 + V 


-20 


1 - Y e 
's 


-20 


(A-4) 
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An effective frequency-dependent input reflection function Y^(jtu) can be 
defined as 


, . % R „ - 
Y 1 J " R s + 


(A- 5) 


Substitution of (A-3) in (A-5) and subsequent simplification yield 


Y i (j«) = 


Y ~ Yr e 
' s L 


-20 


1 ~ VL e 


-20 


(A-6) 


An effective complex output reflection function y 0 (jto) can be defined as 


, v _ R L Z out^ j ^ 
T ° ° " \ + 


(A-7) 


Substitution of (A-4) in (A-7) and subsequent simplification yield 


Y 0 (j<*>) = 


Y, * Y s 
'L 's 


-20 


1 - Vs 5 


-20 


(A-8) 


The quantities y g and y^ relate to reflections of a single wave incident on 
the particular boundry. However, y^Cjuj) and y 0 (jaj) are effective reflection 
functions that take into consideration the effects of multiple reflections. 

Several special cases should be noted. 

(a) Z = R = R 
o o L 


In this case y L = 0, and when this value is substituted in (A-3), (A-6), and 


(A-8), there results 



( A— 10) 


YjQw) = Y g 


V jw) = “V 


-28 


(A-l 1 ) 


Thus, when the image impedance is real and the load end is matched, the input 
impedance is simply the image impedance, and the effective input reflection 
coefficient is simply the source reflection coefficient. 


(b) Z = R = R 
o o s 


In this case y = 0, and when this value is substituted in (A-4), (A-6), and 


(A-8) 


Z _(ja>) = 
out J 


(A-l 1) 


Y 0 ( ja J ) = y l 


(A-l 2) 


Y i (j<*)) = -y L e 


-20 


(A-13) 


Thus, when the image impedance is real and the source end is matched, the output 
impedance is simply the image impedance, and the effective output reflection 
coefficient is simply the load reflection coefficient. 


(O R l = 


In this case, y s = Yl* It can then be shown from (A-6) and (A-8) that 


Y t (jw) = Y 0 C j<*>) = 


- „- 2 ^ V 1 " e 26) 


r L (1 - e " D ) 
. 2 -26 
1 ' V 


, 2 -28 
1 - y e 
's 


(A-14) 
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APPENDIX B 


RECIPROCITY OF TRANSMISSION FUNCTION 

Consider the passive bilateral circuit model shown in figure Bl(a). A source 
with RMS voltage V g and internal resistance R g excites the circuit, and the 
output is connected to a load R^. The output power P Q expressed in terms of the 
RMS load current 1^ is 

p o ■ l l\ 


The maximum available source power P„„ is 

3. V 


P av 4R 

s 

The transmission function |t| 2 is defined as 


(B-2) 



P o _ 'W 

P 2 

av V /4R 
s s 



‘W21 


(B-3) 


where Y 2 ^ is the transfer admittance. 

Assume now that the network is turned around as shown in figure Bl(b). The 
output power P^ is now 

P' = (I t ') 2 R (B-4) 

O L 8 


The available source power is now 


P' 


av 


4R L 


(B-5) 
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The new transmission function 


I t * 1 2 is 


P* 

o 

P' 

av 



4E s R L< Y 21>‘ 


(B-6) 


By reciprocity, however, Y 21 = Y 21 * Comparing (B-6) and (B-3), it can be concluded 
that the transmission function obeys reciprocity for a linear bilateral circuit. 

If, in addition, the network is lossless, the reflection functions | C j cu) | 2 
and |Y 0 (jw)| 2 are expressible as 

I Y i ( j o) ) | 2 = 1 - ! x j 2 = | Y Q ( joi) j 2 (B-7 ) 

since j x | 2 is the same in both directions. Thus, for a linear, bilateral lossless 
network, the reflection functions at both ports are identical. 
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APPENDIX C 


POLES FOR LOW-PASS PROTOTYPE FILTERS 

The following definitions are used in the tables that follow: 

r = pass band ripple in decibels. Table C-l (r = 0) corresponds to the 
Butterworth case, and all other cases correspond to Chebyshev types. 

n = number of poles in low-pass prototype filter. 

= number of resonators in band-pass filter. 

Transmission Constant = numerator of | x( 3 to) | ^ when denominator 
coefficient of s^ n term is unity. 

The value o> = 1 rad/s on a normalized basis is the frequency 
corresponding to the band-edge last ripple bound point. 
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Table C-l. Poles for r = 0 (Butterworth case) 


r = 0 


n 

Poles 

Transmission 

Constant 

2 

-0.7071068 ± jO. 7071068 

1 

3 

-0.5000000 ± jO. 8660 2 54 
-1.0000000 

1 

4 

-0.3826834 ± jO. 9238795 
-0.9238795 ± jO. 3826834 

1 

5 

-0.3090170 ± jO. 9510565 
-0.8090170 ± jO. 5877852 
- 1.0000000 

1 

6 

-0.2588190 ± jO. 9659258 
-0.7071068 ± jO. 7071068 
-0.9659258 ± jO. 2588190 

1 

7 

-0.2225209 ± jO. 9749279 
-0.6234898 ± jO. 7818315 
-0.9009689 ± jO. 4338837 
- 1.0000000 

1 




Table C-2. Poles for r = 0.1 dB 


r = 0.1 dB 


n 

Poles 

Transmission 

Constant 

2 

-1.1861781 ± jl. 3809484 

10.7328410 

3 

-0.4847029 ± jl. 2061553 
-0.9694057 

2.6832106 

4 

-0.2641564 ± jl. 1226098 
-0.6377299 ± jO. 4650002 

0.6708026 

5 

-0.1665337 ± jl. 0803720 
-0.4359908 ± jO. 6677066 
-0.5389143 

0.1677006 

6 

-0.1146934 ± jl. 0565189 
-0.3133481 ± jO. 7734255 
-0.4280415 ± jO. 2830934 

0.04192516 

7 

-0.08384097 ± jl. 0418333 
-0.2349172 ± jO. 8354855 
-0.3394651 ± jO. 4636595 
-0.3767779 

0.01048129 




Table C-3. Poles for r = 0.5 dB 


r = 0.5 dB 


n 

Poles 

Transmission 

Constant 

2 

-0.7128122 ± jl. 0040425 

2.0488703 

3 

-0.3132282 ± jl. 0219275 
-0.6264565 

0.5122176 

4 

-0.1753531 ± jl. 0162529 
-0.4233398 ± jO. 4209457 

0.1280542 

5 

-0.1119629 ± jl. 0115574 
-0.2931227 ± jO. 6251768 
-0.3623196 

0.03201359 

6 

-0.0776501 ± jl. 0084608 
-0.2121440 ± jO. 7382446 
-0.2897940 ± jO. 2702162 

0.008003395 

7 

-0.0570032 ± j 1.0064085 
-0.1597194 ± jO. 8070770 
-0.2308012 ± jO. 4478939 
-0.2561700 

0.002000849 




Table C-4. Poles for r = 1 dB 


r = 1 dB 


n 

Poles 

Transmission 

Constant 

2 

-0.5488672 ± jO. 8951286 

0.9655292 

3 

-0.2470853 ± jO. 9659987 
-0.4941706 

0.2413823 

4 

-0.1395360 ± jO. 9833792 
-0.3368697 ± jO. 4073290 

0.06034558 

5 

-0.0894584 ± jO. 9901071 
-0.2342050 ± jO. 6119198 
-0.2894933 

0.01508638 

6 

-0.0621810 ± iO. 9934115 
-0.1698817 ± jO. 7272275 
-0.2320627 ± jO. 2661837 

0.003771601 

7 

-0.0457089 ± jO. 9952839 
-0.1280736 ± jO. 7981557 
-0.1850717 ± jO. 44294 30 
-0.2054141 

0.0009428962 




Table C-5. Poles for r = 2 dB 


r = 2 dB 


n 

Poles 

Transmission 

Constant 

2 

-0.4019082 ± jO. 8133451 

0.4274285 

3 

-0.1844554 ± jO. 9230771 
-0.3689108 

0.1068571 

4 

-0.1048872 ± jO. 9579530 
-0.2532202 ± jO. 3967971 

0.02671428 

5 

-0.0674610 ± jO. 9734557 
-0.1766151 ± jO. 6016287 
-0.2183083 

0.006678567 

6 

-0.0469732 ± jO. 9817052 
-0.1283332 ± jO. 7186581 
-0.1753064 ± jO. 2630471 

0.001669642 

7 

-0.0345566 ± jO. 9866139 
-0.0968253 ± jO. 7912029 
-0.1399167 ± jO. 4390845 
-0.1552958 

0.0004170907 




Table C-6. Poles for r = 3 dB 


r = 3 dB 


n 

Poles 

Transmission 

Constant 

2 

-0.3224498 ± jO. 7771576 

0.251190075 

3 

-‘0. 1493101 ± jO. 9038144 
-0.2986202 

0.06279750 

4 

-0.0851704 ± .10.9464844 
-0.2056195 ± jO. 3920467 

0.01569938 

5 

-0.0548531 ± jO. 9659238 
-0.1436074 ± jO. 5969738 
-0.1775085 

0.003923656 

6 

-0.0382295 ± jO. 9764060 
-0.1044450 ± jO. 7147788 
-0.1426745 ± jO.2616272 

0.0009812111 

7 

-0.0281456 ± 10.9826957 
-0.0788623 ± 10.7880608 
-0.1139594 ± 10.4373407 
-0.1264854 

0.0002453029 




APPENDIX D 


PROGRAM LISTINGS 


FILTER1 


10 INPUT "NUMBER OF POLES" ;NP 
15 DSHIFT=0 

20 FOR N=1 TO NP STEP 2 

30 INPUT "TRANS DENOM: SIGMA, OMEGA(+VALUES)" ; SIGMAD(N) ,OMEGAD(N) 
40 SIGMAD(N)=-SIGMAD(N) 

50 NEXT N 
55 T2MULT= 1 

60 FOR N=2 TO NP STEP 2 
70 SIGMAD(N)=SIGMAD(N-1) 

8 0 OMEGAD ( N ) =-OMEGAD ( N- 1 ) 

90 NEXT N 

100 INPUT "TRANS CONSTANT" ; T2CON 

110 INPUT "VALUE OF D";D 

120 INPUT "FREQ INCREM" ; DELW 

130 INPUT "INTEGER OF HIGHEST FREQ";KMAX 

132 FOR N=1 TO NP 

134 SIGMAD(N)=SIGMAD(N)+DSHIFT 

136 NEXT N 

138 T2CON=T2CON*T2MULT 
140 LPRINT "POLES" 

150 FOR N=1 TO NP 

160 LPRINT SIGMAD(N) ,OMEGAD(N) 

170 NEXT N 
180 LPRINT 

190 LPRINT "NUMBER OF POLES =";NP,"D =";D,"NUM CON =”; T2CON 
200 LPRINT 

210 LPRINT "FREQUENCY" TRANS SQ","ACCUM TRANS" ,"ACCUM 4" 

220 LPRINT, "DB RESPONSE" , "ACCUM NOR TRAN" , " STAT NOR BW" , "DEGRAD" 

230 W=0 

240 INTTSQ=0 

250 INTT4=0 

260 A2=0 

270 A4=0 

280 TRANS=0 

290 BSTATN=0 

300 DEGRAD=0 

310 FOR K=0 TO KMAX 

320 DF1=1 

330 FOR N=1 TO NP 

340 DF2=(SIGMAD(N)-D)*2+(W-OMEGAD(N)) A 2 

350 DF1=DF1*DF2 

360 NEXT N 

370 TSQ=T2C0N/DF1 

380 TDB=4.34294*LOG(TSQ) 

390 T4=TSQ A 2 

400 IF K=0 THEN 460 

4 1 0 INTTSQ= INTTSQ+ . 5 * ( TSQ+A2 ) * DELW 
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420 INTT4=INTT4+ . 5* ( T4+A4 ) *DELW 

430 TRANS= INTTSQ/ W 

440 BSTATN=INTTSQ A 2/( INTT4*W) 

450 DEGRAD=1/(TRANS*SQR(BSTATN)) 

460 LPRINT W,TSQ, INTTSQ, INTT4 

470 LPRINT, TDB, TRANS, BSTATN, DEGRAD 

480 W=W+DELW 

490 A2=TSQ 

500 A4=T4 

510 NEXT K 

520 LPRINT 

530 LPRINT 

540 INPUT "ANOTHER RUN (Y OR N)";B$ 

550 IF B$=" N" THEN 630 

551 FOR N=1 TO NP 

552 SIGMAD(N)=SIGMAD(N)-DSHIFT 

553 NEXT N 

554 T2C0N=T2C0N/T2MULT 

555 DSHIFT=0 

556 T2MULT= 1 

560 INPUT "DO YOU WANT TO PREDISTORT (Y OR N)";C$ 

570 IF C$="N" THEN 110 

580 INPUT "VALUE OF D SHIFT" ; DSHIFT 

590 C1=SIGMAD(1) A 2 

600 C2=(SIGMAD( 1 )+DSHIFT) A 2 

610 T2MULT=C2/C1 

620 GOTO 110 

630 END 



FILTER2 


10 INPUT "NUMBER OF POLES" ;NP 
1 5 DSHIFT=0 

20 FOR N=1 TO NP STEP 2 

30 INPUT "TRANS DENOM : SIGMA, OMEGA( +VALUES ) " ; S IGMAD(N ) , OMEGAD( N ) 
40 SIGMAD(N)=-SIGMAD(N) 

50 NEXT N 
55 T2MULT=1 

60 FOR N=2 TO NP STEP 2 
70 SIGMAD(N)=SIGMAD(N-1 ) 

80 OMEGAD ( N ) =-0MEGAD( N- 1 ) 

90 NEXT N 

100 INPUT "TRANS CONSTANT" ;T2CON 

110 INPUT "VALUE OF Q";Q 

115 D=l/Q 

120 DELW=.0 1 

132 FOR N=1 TO NP 

134 SIGMAD(N)=SIGMAD(N)+DSHIFT 

136 NEXT N 

138 T2CON=T2CON*T2MULT 
140 LPRINT "POLES" 

150 FOR N=1 TO NP 

160 LPRINT SIGMAD(N) ,OMEGAD(N) 

170 NEXT N 
180 LPRINT 

190 LPRINT "NUMBER OF POLES =" ;NP, "Q=" ;Q, "NUM CON =";T2CON 
200 LPRINT 

210 LPRINT "FREQUENCY", "MODIFIED Q, "DEGRAD" 

220 RCOMP=-10 
230 W=0 
235 TDB=0 
240 INTTSQ=0 
250 INTT4=0 
260 A2=0 
270 A4=0 
275 ADEG=0 
280 TRANS=0 
290 BSTATN=0 
300 DEGRAD=0 
310 FOR K=1 TO 2 
320 DF1=1 
330 FOR N=1 TO NP 

340 DF2=(SIGMAD(N)-D) A 2+(W-OMEGAD(N) ) *2 

350 DF1=DF1*DF2 

360 NEXT N 

370 TSQ=T2C0N/DF1 

380 TDB=4. 34294 *LOG(TSQ) 

390 T4=TSQ A 2 

4 10 INTTSQ=INTTSQ+ . 5* (TSQ+A2 ) *DELW 
420 INTT4=INTT4+.5*(T4+A4)*DELW 
430 TRANS=INTTSQ/W 
440 BSTATN=INTTSQ A 2/(INTT4*W) 

450 DEGRAD=1/ (TRANS*SQR(BSTATN) ) 

460 IF TDB<=RCOMP THEN 504 
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470 ADB=TDB 
480 W=W+DELW 
490 A2=TSQ 

500 A4=T4 

501 ADEG=DEGRAD 

502 GOTO 320 

504 ADJ=(RCOMP-ADB) / (TDB-ADB) 

506 WMOD=W-DELW+ADJ*DELW 

508 DMOD=D/WMOD 

509 QMOD=l/DMOD 

510 LPRINT WMOD,QMOD, DEGRAD 
512 RCOMP=RC0MP-10 

514 ADB=TDB 
516 W=W+DELW 
518 A2=TSQ 

520 A4=T4 

521 ADEG=DEGRAD 

522 NEXT K 
525 LPRINT 
530 LPRINT 

540 INPUT "ANOTHER RUN (Y OR N)";B$ 

550 IF B$="N" THEN 630 

551 FOR N=1 TO NP 

552 SIGMAD(N)=SIGMAD(N)-DSHIFT 

553 NEXT N 

554 T2CON=T2CON/T2MULT 

555 DSHIFT=0 

556 T2MULT=1 

560 INPUT "DO YOU WANT TO PREDISTORT (Y OR N)";C$ 

570 IF C$="N" THEN 110 

580 INPUT "VALUE OF D SHIFT" ;DSHIFT 

590 C1=SIGMAD(1) A 2 

600 C2=(SIGMAD( 1 )+DSHIFT) A 2 

610 T2MULT=C2/C1 

620 GOTO 110 

630 END 



APPENDIX E 


CURVES OF DEGRADATION FACTOR VS BAND-EDGE ATTENUATION WITH q AS A PARAMETER 

All curves contained in Appendix E are plots of the degradation factor vs band- 
edge attenuation. To simplify the presentation, each caption provides only the pass- 
band ripple and the number of poles for the particular set of curves. 

In general, the following low-pass reference q values were considered: », 

100, 30, and 10. The values of q were established at the lossless low-pass proto- 
type "cutoff" frequency corresponding to the filter type. For Butterworth filters 
(r = 0 dB) and 3-dB Chebyshev filters (r = 3 dB), this corresponds to the frequency 
at which the response is down 3 dB. For other ripple levels, it is the highest 
frequency at which the response crosses the ripple bound. 

For the three lossy q values, both standard and predistorted designs were con- 
sidered. For lower order filters and higher values of q, the standard and pre- 
distorted cases essentially coincided. For higher order filters and lower values of 
q, the curves were off the scales and are not shown. 

The low-pass reference q values are denoted in a row above, and in the same 
sense as, the curves. 

The pole values in Appendix E were not adjusted to compensate for shifts in 
reference frequencies as the dissipation increased. 

NOTE : In all cases, solid curves represent standard designs, and broken lines 

represent predistorted designs. 
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APPENDIX F 


CURVES OF DEGRADATION FACTOR VS q WITH RIPPLE AS A PARAMETER 
All curves contained in Appendix F are plots of the degradation factor vs low- 
pass reference q for two specific values of band-edge attenuation: 10 dB and 

20 dB. To simplify the presentation, each caption provides only the attenuation and 
the number of poles for the particular set of curves. 

The poles were adjusted by the computer program so that the q values obtained 
in the data represent the actual low-pass q at the frequency at which the attenua- 
tion is the specific level (i.e., 10 dB or 20 dB). 

NOTE: In all cases, solid curves represent standard designs, and broken lines 

represent predistorted designs. 
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APPENDIX G 


SUMMARY CURVES 

A summary of the minimum degradation factor that can be obtained (with the 
filter types investigated) as a function of the number of poles for 10 dB specified 
bandedge attenuation is provided in figure G1 . A similar summary for 20 dB 
specified bandedge attenuation is provided in figure G2. 

The range of minimum degradation factor based on the bounds of three resonators 
and seven resonators is provided in figure G3 . 
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